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Abstract 

The non-Abelian analog of the classical Coulomb gas is discussed. The statistical 
mechanics of arrays of classical particles which transform under various representa- 
tions of a non-Abelian gauge group and which interact through non-Abelian electric 
fields are considered. The problem is formulated on the lattice and, for the case of 
adjoint charges, it is solved in the large N limit. The explicit solution exhibits a first 
order confinement-de-confinement phase transition with computable properties. In 
one dimension, the solution has a continuum limit which describes 1+1-dimensional 
quantum chromodynamics (QCD) with heavy adjoint matter. 
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1 Introduction 



The classical Coulomb gas is an important model in statistical mechanics. It is solvable 
in one dimension and in two dimensions it exhibits interesting critical phenomena. In this 
Paper, we shall formulate a non-Abelian generalization of the Coulomb gas. We consider 
the thermodynamic properties of arrays of non-dynamical quarks which transform under 
irreducible representations of the gauge group and which interact via non-Abelian electric 
fields. The main motivation is to study the confinement- deconfinement phase transition 
which is expected to occur in gauge theory at sufficiently high temperature or density. 

The model which we formulate is most interesting in the case where the quarks trans- 
form under the adjoint representation of the gauge group. In that case, it can be solved 
explicitly in the large N limit and exhibits a first order phase transition when the spa- 
tial dimension is D=l,2,3. When D=l, it has a continuum limit where it represents 
(1+1) -dimensional QCD with very heavy adjoint matter fields. 

In the absence of quarks, i.e. the zero density limit, the model which we consider 
reduces to the strong coupling limit of lattice Yang-Mills theory. In that limit, only the 
color-electric interactions are retained. The phase structure of that system was studied in 
the large N limit in refs.[[I], 0. There, it was found that in D=2 and D=3, there is a first 
order phase transition which occurs as the parameter 7 = e 2 N/2T is varied. The strong 
coupling, low temperature, confining phase occurs at 7 >> 1 and the weak coupling, high 
temperature, de-confined phase at 7 << 1. In the former case one can find an explicit 
solution of the lattice theory in the large N limit. This phase is stable to small fluctuations 
for 7 greater than some critical value. In the small 7 limit, approximate techniques are 
required and a de-confined phase is found. There is a co-existence region where both 
phases are stable to small fluctuations, and are therefore separated by an energy barrier. 
In that region, there is a first order phase transition between them. 

In D = 2, 3, the strong coupling theory is not renormalizable. The phase transition is 
first order, and there is no limit where the latent heat is small, i.e. of order N 2 times a 
quantity which remains finite as the lattice spacing is taken to zero. One could speculate 
that, if the strong coupling limit were corrected to include the magnetic term in the 
Hamiltonian, one would recover a renormalizable theory, manifest in the fact that there 
would be some limit in which the first order phase transition had finite latent heat in the 
continuum. There are good arguments for this scenario in the literature. As well, large 
N techniques have been used to study the confinement-deconfinement transition in pure 
Yang-Mills theory @. 

In D=l, there is no magnetic term in the full Hamiltonian, and Yang-Mills theory can 
be solved exactly using a technique similar to the strong coupling approximation. Rather 
than being generated dynamically, the confining quark-antiquark potential appears at the 
tree level where the string tension is proportional to the coupling constant e 2 . However, 
in D=l, pure Yang-Mills theory is trivial in that it has no propagating degrees of freedom. 
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It therefore does not exhibit a phase transition when the temperature is non-zero and it 
is in the confining phase for all values of 7.^ 

In this paper, we consider the model where a gas of charged sources (quarks) in the 
adjoint representation of the gauge group are added to strongly coupled Yang-Mills theory. 
For D = 2, 3 there is still a first order phase transition which occurs on a critical line in 
the 7 — A plane, where A ~ e _At//T is the fugacity and /1 is the chemical potential of the 
quarks. When 7 is greater than a certain critical value, the deconfinement transition can 
be induced by increasing the density, which is controlled by increasing A. There is also 
a third order phase transition in the de-confined phase, similar to the Gross- Witten || 
transition of D=l lattice Yang-Mills theory. 

In D=l, the coupling of sources to Yang-Mills theory makes the model non-trivial. 
The continuum version of this model has been studied in ref.0. There is a first order 
phase transition between the confining and de-confined phases which can be obtained by 
increasing A for any value of 7. This phase transition originates in a percolation transition 
for electric flux lines. The quantum states are arrays of adjoint quarks on the line with 
non-dynamical strings of electric flux joining them. Each quark must have one line of flux 
entering and one line leaving it. For a fixed number of quarks, the model is explicitly 
solvable and has a finite dimensional Hilbert space, corresponding with the different ways 
of distributing electric flux so that the states are gauge invariant. In the large N limit, 
the energy of a state is proportional to the total length of lines of electric flux plus the 
chemical potential times the number of quarks. At low temperature and density, the 
statistical sum is dominated by configurations which are a dilute gas of mesons - color 
neutral bound states of two or more adjoint quarks bound together by the appropriate 
number of strings of electric flux. 

The property of confinement is characterized by measuring the energy necessary to 
insert a fundamental representation quark- ant iquark pair into the system. Gauge in- 
variance requires that the pair is connected by a single string of electric flux. At low 
temperature, the energy of this string is proportional to its length, leading to the linear, 
confining interaction. In a deconfined phase, for large separation, the energy would be a 
constant as the distance between the quark and antiquark is varied. This occurs when, at 
the deconfinement phase transition, the strings percolate in the one-dimensional space, 
and the addition of a flux string between a quark and an antiquark adds a small amount 
(or no) energy to the energy of the typical configuration. 

In the next few subsections, we review some of the formalism which we require in the 
later sections - the Hamiltonian formalism of lattice Yang- Mills theory, the construction of 
the partition function in that formalism and some properties of states in the non-Abelian 
Coulomb gas. 

3 1+1 dimensional Yang-Mills theory on the sphere is known to have a third order Kazakov-DouglasQ 
phase transition. However, in the cylindrical geometry which is appropriate to finite temperature Yang- 
Mills theory on the open line, R , there is no third order phase transition. 
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1.1 Hamiltonian formulation of lattice Yang-Mills theory with 
external sources 



We begin by reviewing the Hamiltonian formalism of Yang-Mills theory on the lattice. 
This is standard material which can be found in ref.0, for example. In the Hamiltonian 
formalism, time is continuous and the space is approximated by a hypercubic lattice with 
sites x, y, . . . and oriented links I. The spatial gauge fields are operator valued unitary 
matrices, Ui, and electric fields are operator valued Lie algebra elements, Ei, 



U t U} = 1 = Ufa x , E} = E l 
Both are associated with oriented links and obey the reflection conditions 



U-i = U\ , E^ 



-UjEiUi 



The electric fields can be expanded in a basis 

E l = Y,EfT A 



1.2) 



;i.3) 



where T A are hermitean generators in the fundamental representation of the Lie algebra 



T , T 



B 



if 



ABCrpC 



The operators have the algebra 



[17,, 17,,] = 
The Hamiltonian is 



E A ,E*]=if ABC E?5 lv 



E?,Uv 



T USi, 



^=EyK) 2 +E^Tr(n^+n^) 

LA □ Ze V □ □ / 



;i-4) 



;i.5) 



;i.6) 



and the gauge constraint is 



K 



ldn(x) 



(1.7) 



i=l 



The first term in the Hamiltonian is the electric energy and the second term, which is 
summed over oriented elementary plaquettes □, is the magnetic energy. In Q (x), n(x) is 
the set of links one of whose endpoints is the site x and with orientation toward x. Also, 
we have considered the system in the presence of an array of K classical quarks which 
transform under representations Ri and are situated at lattice sites x%. 

The gauge constraint commutes with the Hamiltonian and generates the gauge trans- 
formation, 

U x - Uf = g{x)U l9 \y) , E t - Ef = g{x)E l9 \x) (1.8) 
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where g(x) is a unitary matrix in the fundamental representation of the gauge group and 
x and y are the endpoints of the oriented link I: 51 = [y] — [x]. It is useful to consider the 
Schrodinger representation where quantum states are functions of the gauge fields and the 
action of the electric field on these states is defined algebraically. In this representation, 
the gauge constraint implies that physical states transform as 

^...a K p] = g^M) ■ ■■9aXM ip bl ... bK [u 9 } (1.9) 

Since the Hamiltonian is gauge invariant, the eigenstates of the Hamiltonian carry a 
representation of the gauge group. The physical states are those which transform like 

In the continuum limit, the spatial gauge field A is obtained from U\ = Ve 1 ^ A ' dl ~ 
1 + ial • A + . . . where a is the lattice spacing and I is a unit vector in the direction of 
the link I. Also, the continuum electric field operator is defined as Ef = a D_1 / • E A . 
In D > 1, the magnetic energy reduces to the magnetic field squared. The Hamiltonian 
(|1.6|), gauge constraint ( |1.7| ) and operator algebra ( |1.5| ) reduce to those of continuum 
Yang-Mills theory, 

H = jd D x(^{E A Y + ^{B A )^ (1.10) 
G A {x) = (D.E) A + j: Ti8{x- Xi ) ~ (1.11) 

i=l 

and 

[A, l (x),A 1 (y)] = , [E i (x),E j (y)] = 

respectively. On the lattice, energies are measured in units of the lattice spacing. It is 
only the very low-lying states of the dimensionless lattice Hamiltonian, with energy of 
order a ■ e which have finite energy, e, in the continuum limit, a — > 0. 



Ef(x),Af(y) =-i6 A "6 ij 6(x-y) (1.12) 



AB • 



1.2 Partition function 

The partition function at temperature T is the trace of the Gibbs distribution e~ H l T over 
gauge invariant physical states. It is convenient to take this trace in the "coordinate" 
representation. The complete set of states is obtained by taking the direct product of the 
"position" eigenstates of the unitary matrices 

\U)=U\Ui) (1-13) 
i 

for each link I, together with some basis elements, e ai . . . e aK , which carry the representa- 
tions Ri . . . Rk of the gauge group. The states \Ui) are normalized by 

(U l \U' l ) = 5(U l ,U' l ), (1.14) 
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where 5(Ui,U{) is invariant ^-function on the group manifold. The states (|1.13|) do not 
satisfy Gauss law constraint (|1.7|). The projection on the physical subspace can be realized 
by gauge transforming the states at one side of the trace and subsequently integrating 
over all gauge transforms. The result is 

Z[x l ,R i ,T] = f Yl[dUi\Y[[dg(x)] (U\e- H ' T \UV) Trg*( Xl ) . . . Trg R «(x K ) (1.15) 

J l x 

where [dUj\ and [dg(x)] are invariant Haar measures. We consider the case where all 
external quarks are in the adjoint representation. In the adjoint representation, the trace 
can be taken as the modulus squared trace of the fundamental representation group 
element, 

Tr# adj (x) = \Tig{x)\ 2 -l (1.16) 

We then multiply by the K'th power of the fugacity, \ K , multiply by l/K\, sum over 
positions, Xi, . . . , xk and sum over K. This produces the effective theory 

Z[\T) = f UldUt] U[dg(x))e- S °«™ (1.17) 

J l X 

where the effective action is 

e -S es[ U,9] = {Ule -H/T lug)e J2 x X(\Tr g (x)f-l) (Llg) 

This effective action has gauge symmetry, 

S eS [hUh\hgtf] = S cS [U,g] (1.19) 

Where h is an element of the gauge group. It also has a global symmetry 

S eS [U,zg] = S cS [U,g] (1.20) 

where z is a constant element of the center of the gauge group. This latter symmetry is 
related to confinement - |T2[ . For the gauge group SU(N), the center is and it is 
referred to as Zjv-symmetry. For U(N), the center is U(l). 

When this symmetry is realized faithfully, the theory is confining. When it is spon- 
taneously broken, the system is in the de-confined phase. An order parameter for this 
symmetry is the Polyakov loop operator^ which is the trace of the group element, 

P(x) = Trg(x) (1.21) 



4 In the spacetime path integral formulation of finite temperature gauge theory |13| Aq(t, x) is a 
Lagrange multiplier field which enforces the gauge constraint. The Euclidean action is 



J Q 1/T dT (e EiU i^ Ui - H ^ E]+ % A ° (r ' x)gA{T ' ■ 

Z = J dA Q dE[dU]e 



S = J dr | > E,r'<—I' i -!l\(-.E}-\- l y A;}[-..r)G A iT..r) 
and the partition function is 



-S[A ,E,U] 
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where g(x) is in the fundamental representation of the gauge group. It is gauge invariant 
and, under ( |1.20| ), it transforms as 



P(x) -> z P(x) (1.22) 

The quantity 

F(x,R,\,T) = -T\n(P(x)} (1.23) 

is the free energy which is necessary to insert a classical source in the fundamental rep- 
resentation into the system. In the confining phase, the expectation value vanishes and 
this free energy is infinite. In a de-confined phase the expectation value is non-zero and 
the free energy is finite. 

The two-point correlators of Polyakov loops are related to the free energy of the 
quark-ant iquark pair inserted in the vacuum. Assuming that the cluster property holds, 
the two-point correlator behaves at \x — y\ — > oo as 

(P(x) Pt(y)) — > | (P{x))\ 2 + _ C 7 °^_ 1)/2 e -*\-* + ..., (1.24) 

where M\ is the mass of the lowest excitation with appropriate quantum numbers. The 
two-point correlator goes to zero if the Polyakov loop expectation value vanishes and goes 
to a constant otherwise. In the deconfining phase, the quark-antiquark potential can be 
defined by subtracting the self-energies (|1 . 23| ) from the logarithm of the correlator ( |1.24 ): 

V(x,y) = -Tln((P(x) P j (y))) +T\n(P(x)} +Tln(P t (y)) 

( ' OUSt .e-"*l-*l. (1.25) 



\ x _j,|(-D-l)/2 

In the confining phase the self-energies are infinite, but the potential still can be defined. 
It grows linearly with distance: 

V(x,y) = -Tin ((P(x) P^(y))) ~ a\x - y\, (1.26) 

where the string tension is 

a = TMi. (1.27) 

For the models that we shall consider, the mass M\ will be calculated exactly in Sec. |2] and [| 

with periodic boundary conditions in Euclidean time. In the strong coupling limit, this model can be 
reduced to ( 1.15 ) by integrating the time dependent modes explicitly. The dynamical variable g(x) is 
the holonomy group element for transport of a fundamental representation charge around the periodic 
Euclidean time, 

g(x) = Pexp^i J drA (r,x)^ . 
The Polyakov loop operator is the trace of this holonomy element P(x) = Trg(x). 
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1.3 Adjoint non-Abelian Coulomb gas 



We shall treat the effective action (JTTH) exactly in D = 1 where there is no magnetic 
term in the Hamiltonian. In D > 1, we consider the strong coupling limit where 

e 2 ^oo , T^oo , 7 = e 2 N/2T finite (1.28) 

(We shall eventually also consider the limit where iV — > oo with 7 finite.) Here e 2 and T 
(as well as the lattice Hamiltonian H) are dimensionless quantities which will eventually 
get their engineering dimensions from factors of the lattice spacing. In this limit, the 
magnetic term in the Hamiltonian is ignored, so that the Hamiltonian is given by the 



electric term only and the Boltzmann weight (|1.18| ) factorizes on the product of the heat 
kernels 

K(g,h\r) = (h\e^ A / N \g) (1.29) 

for the Laplacian 

A = £(V) 2 (1.30) 

A 

on the group manifold with "time" 7 = e 2 N/2T. In this limit, the magnetic interactions 
of the external quarks as well as the magnetic self-interactions of gluons are absent. It 
isolates the electric interactions and self-interactions. 

In this strong coupling limit, the partition function has the form 

z = /ntwin^i n #(^,^17) e ^^ 2 ^ (im) 

J x I I 

It is this model which we call the adjoint non-Abelian Coulomb gas. 

1.3.1 Lattice string statistical mechanics 

In strong coupling limit, the effective Hamiltonian is 

5» = Ei(^ (1-32) 

l,A / 

The local moments of the electric field distribution are conserved, 



E A (x),H 



(1.33) 



and distribution of electric fields is frozen in time. The eigenstates of the Hamiltonian 
( |1.32j ) can be constructed by acting with link variables Ui on the vacuum. On a given 



link, we begin with the singlet state, 

Ef\0) = (1.34) 
and the action of Ef on excited states is defined combinatorially, 

Ef (UiU |0> = (T^U |0> (1.35) 
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Figure 1: Graphical representation of the states \^f±) 

Ef (UiUiUi)^ |0> = {T A Ui) ab (Ui) cd |0) + (U^T^U |0) (1.36) 

and so on. The gauge invariant states contain closed loops of electric flux and lines of 
electric flux, appropriate numbers of which end at the quarks. Closed lines of flux are 
created by the traces of the spatial Wilson loop operators in different representations of 
the gauge group, for the closed curve T, 

W R [r\=TrHU l R (1.37) 

For a non-intersecting array of electric flux strings, the energy is proportional to the total 
length of all strings, 

E = ^C 2 ^L[T] (1.38) 
z r 

Strings interact when they have common links. For example, we consider an array of 
strings where two intersecting strings share a common link. In this case, the eigenstates 
of the Hamiltonian are linear combinations of the uncrossed and crossed strings, fig. [If 

|*±) = {UabUcd ± U ad U cb ) |0> . (1.39) 

The contribution of the states \^±) to the energy of the string configuration is equal to 
e 2 {N - 1)(JV + 2)/N and e 2 {N + 1)(JV - 2)/N, respectively. In the large N limit, these 
states are degenerate and diagonalization of the Hamiltonian does not mix the crossed 
and uncrossed strings. Therefore, the interaction vanishes in the large N limit. Thus, in 
the large N limit, the statistical mechanics of the strong coupling Yang-Mills theory is 
equivalent to a statistical model of non-interacting lattice strings for which the partition 
sum can be written 

z = Y J e~ c ^ mi2N = J2 n ( L ) e ~ C2lL/2N - ( L4 °) 

r l 

Such a system is expected to have a percolation transition when the entropy of configu- 
rations of strings overtakes the energy fI3}| . For large curves, the number of strings with 
a given length L grows like 

n(L) ~ const. (2D - 1) L (1.41) 

The critical temperature is given by 

7crit . = 21n(2D-l) (1.42) 

This is the result which was found in [|], 3 by the direct solution of the model. 

When there are adjoint quarks present, quantum states still contain closed loops of 
electric flux. In addition, it is necessary that each quark absorbs and emits one line of flux. 
Increasing the density of quarks decreases the phase transition temperature somewhat. 
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1.3.2 Continuum limit 



In the continuum limit, the trace which is taken to obtain the partition function can be 
taken over a complete set of eigenstates of the gauge field operator, twisted by the gauge 
transformation [15|], (A \ e~ H ^ T \A 9 ), where A 9 = g(A + id)g^ . Formally, we can consider 



the strong coupling limit and retain only the electric term in the Hamiltonian: 

(A\ e -^S {EA)2 \A 9 ) = const, exp (- J J d D xTr(A - A 9 ) 2 ^ 

= const, exp (- — J d D xTr \Dg(x)\ 2 ^j (1.43) 
where Dg = Vg + i [A, g]. The partition function is 

Z = J[dA][dg] e -/^-[^I^WI 2 -A(l^)l 2 -D] (L44) 

The effective action is that of a gauged principal chiral model in D dimensions. The con- 
tinuum treatment of the strong coupling limit, however, is valid only at D = 1, where the 
magnetic term in the Hamiltonian is absent from the beginning; this model is considered 
in the next section. At D > 1 the field theory defined by ( |1.44| ) is nonrenormalizable 
unless the kinetic term for the gauge fields is added. 



The mean density of quarks is obtained from the partition function as 

Y- = ^A^lnZ[A,T]=A(Tr^-l)^ A \(Trg) | 2 + 0(A°) (as N - oo) (1.45) 

where we have assumed that the expectation value is translation invariant. The last 
equality follows from factorization of invariant correlators in the large N limit. From eq. 
( P-.45Q we see that, in the large N limit, the Coulomb gas in the confining phase where 
| (Tr g(x)) | = | (P(x)) | = is dilute with density is of order one. This is a result of the fact 
that the number of species of glue-balls and hadrons does not grow in the large N limit. 
Above the deconfining temperature, where \{P(x))\ ^ there are N 2 — 1 gluon and N 2 
quark degrees of freedom and the density is of order N 2 . A similar consideration applies 
to the free energy. The free energy in the confining phase should grow like the number of 
degrees of freedom and be O(N ) in the large N limit. In the de-confined phase the free 
energy should be of order N 2 . For N infinite, this implies that a phase transition between 
the two phases would have infinite latent heat. However, one should consider the large 
iV limit as analogous to the infinite volume limit in statistical mechanics where, strictly 
speaking, there is no phase transition in finite volume, but the analysis in infinite volume 
is to a good approximation accurate in the physical finite system too. 

This will be particularly relevant to the case of D = 1, where, if N is finite, the 
dimensionality of the system is too low to have spontaneous symmetry breaking. It has 



10 



local interactions and is effectively a one-dimensional system. This implies that, when N 
is finite, it is always in the confining phase with unbroken center symmetry. In the infinite 
N limit, there can be a phase transition and a phase with broken symmetry. However, 
the latter limit should describe the physical behavior of the system accurately even when 
N is finite if the "finite size" corrections of order 1/N 2 are small. 

Even this case is subtle, since the effective symmetry in the infinite N limit is U(l), 
a continuous symmetry and the effective dimension is two. The small 7 phase can be 
ordered only because the action is non-local in index space, similar to an infinite-ranged 
spin model. 

In Section 2 we shall analyze the large N limit of the one dimensional model directly 
in the continuum. We find an exact solution of the continuum model and show explicitly 
that it has a first order confinement-deconfinement phase transition which occurs at a 
critical line in the 7 — A plane. In Section 3 we discuss the strong coupling lattice model 
in any dimensions. We find an exact solution in the confined phase and an approximate 
solution of the deconfined phase. As on one dimension, there is a first order confinement- 
deconfinement transition. In Section 4, we discuss the results. 



2 Large N limit of one— dimensional model 

It is possible to analyze the large N limit of the D=l model directly in the continuum 
limit. We consider the gauge group U(N) which has center £7(1), or SU(N) with center 
Zn - there is no difference in the large N limit. In the gauged principal chiral model 
( |1.44j ), since the space is an open line, we can choose the gauge A = 0. The partition 
function for the resulting model, 




is equivalent to that of unitary matrix quantum mechanics (where x is imaginary time) 
and can be solved in the large N limit by the methods of collective field theory [[II], |T^| . 
The method is essentially based on the relation between matrix quantum mechanics and 
nonrelativistic fermions ||18|| . 

The dynamical variables in (|2.1|) are the phases ajt(x) of the eigenvalues e™ 1 ' 1 ' of the 
Polyakov loop variables g(x). They can be interpreted as coordinates of fermions which 
live on a circle. In the large N limit, the statistical integral in (|2.1|) is dominated by a 
single distribution of eigenvalues. We introduce the eigenvalue density 

1 N 

p{e,x) = -Y,^-a k {x)) 1 (2.2) 

k=l 

which is a periodic function which is normalized to unity on the interval (—71,71). It has 
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the mode expansion 



p(6,x) = ±-\l + J2c n (x)e in0 ] , < = c_ n (2.3) 



In a truly confining phase, all moments of the eigenvalue distribution vanish, p con f. = 1/27T. 
If any of the other moments are non-zero, there are Polyakov loop operators which have 
non-zero expectation values, if c n (x) ^ then (Tr g n (x) ) 7^ 0. Thus, all combinations of 
c n ( x ) 7^ characterize all de-confined phases. 

In the large N limit the eigenvalue density obeys a classical, saddle-point equation 
which can be deduced from canonical analysis of the collective field theory Hamiltonian [|T6f , 
0: 



H= dO 



A 



d9p(9)e 



2 7 

- — (2-4) 
24 v ; 



with subsequent Wick rotation to an imaginary time. Here H(9, x) is the variable which 
is the canonical conjugate to p(9,x), so that the Poisson bracket is 

{ P (9,x),Il(9',x)} = 6(9-9') (2.5) 

(here x is the time variable) and v{9) = dU/d6 is the velocity of the Fermi fluid. The 
equations of motion, following from ( |2.4| ), read, after the change x — > ix, v — > —iv, as 
follows: 

| + 7^) = 0, (2.6) 

I + V>% ~ - 2 1P% + 2AIm = 0, (2.7) 

where 

Cl ( x ) = f d6p(6,x)e id . (2.8) 



It is expected that the solution of these equations corresponding to the physical vacuum 
is a constant po{9). 

One may expect that, at least at sufficiently low temperature or, equivalently, at 
sufficiently large 7, the system is in the confining phase with unbroken center group 
symmetry. This symmetry acts on p(6, x) by translation 9 — > 9 + 9q, so the only symmetric 
solution is p C onf.(#) = It always satisfies the equations of motion, as Ci(x), defined by 
eq. ( |2.8| ), is equal to zero when p = p con f.- 

However, this solution, which is an exact solution of (|2.6| ) and (|2.7| ), is stable to small 
fluctuations only if 7 is large enough. It becomes unstable for 7 < 7 C (A). To see this, it is 
necessary to analyze the spectrum of fluctuations in the strong coupling phase. Consider 
p(9,x) = 2^(1 + ip(9, x)) and consider the equations ( |2.6p and ( |2.7] ) to linear order in (p: 
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lH! + 2AIm(ci e-) = 0. 
Differentiating eq. ( |2.9| ) with respect to x and eq. ( |2.10| ) with respect to 9 and subtract- 
ing the latter, multiplyed by 7, from the former, we obtain an equation for the density 
fluctuations: 

The eigenmodes are the Fourier harmonics of <p(9,x): 

v(9,x) = Y,Cn(x)e™ e (2.12) 



where now c n are infinitesimal: 



— A7(5 n>1 + 8 n -i) 



n^O. (2.13) 



We obtain the following spectrum of excitations: 



MZ = l^-\>y8 nl , n=l,2,... (2.14) 

At 7 = 7c (A), where 

7c (A)=4A, (2.15) 

the lowest eigenvalue (n = 1) goes to zero. For smaller 7 this eigenvalue is negative 
and leads to the instability of the strong coupling solution where c\ is the first mode 
to become unstable. However, for reasons which will become clear once we consider the 
weak coupling phase, 7 C (A) should not be identified with the point of the deconfining 
phase transition. 

The contribution of fluctuations to the free energy is given by 

SF/V ^r* fe> (e + ff) + ,„ (q^i)) (2,e) 



00 



2vr \V V 4 J \ k 2 + 7 2 /4 



the singular part of which, near the critical line is SF/V ~ ^f- (Jl-AX/^ - l). 

The deconfining solution can be obtained by integration of eq. (|2.7|) at v = 0. The 
density Po(0) can always be chosen to be an even function of 9. Thus c\ is real, and one 
finds from eq. (|2.7|) : 

p (9) = -J-(e + 2\ Cl cos 9). (2.17) 



7T V 7 

The Fermi energy E and the constant c\ are to be determined from the normalization 
condition and eq. 



/•tfmax d9 I 2 / \ 

/ — 4 /-(£? + 2Aci cos 0) = 1, (2.18) 
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A 



Figure 2: The large N phase diagram of the one-dimensional model. I - strong coupling 
(confining) phase, II - weak coupling (deconfining) phase; 1 - line on which the weak 
coupling phase terminates: 7 = 7* (A), 2 - line of the first-order phase transition: 7 = 
70(A), 3 - line of the instability of the strong coupling phase: 7 = 7 C (A) 



/ — cos0W-(£ + 2AciCos0j = c h (2.19) 



M „ 12 

7 

6*max = 7r - arccos — — . (2.20) 

It follows from these equations that fl mffl tends to zero at 7 — ► and grows with the 
increase of 7. Eventually it reaches 7r, where the weak coupling phase terminates, because 
the eigenvalue distribution begins to overlap with itself due to 27r-periodicity. At the 
critical point E* = 2Aci*, the integrals in (|2.18|) and Q2.19|) can be done explicitly, and we 
find that c u = ~ and 

7 ,(A) = ^A = 4.323A. (2.21) 

At this critical line, the weak coupling phase is unstable. It is there that the translation 
zero mode of the linearization of the equations (|2 . 6| ) -def c , dp/ 89 becomes normalizable and 
the resulting fluctuations restore the translation invariance in 9, i.e. the center symmetry. 
In the phase with density given by wcs this zero mode is not normalizable and therefore 
it is ineffective in restoring the symmetry. 

We obtain the following picture of the deconfining phase transition (fig. Q). The weak 
and strong coupling phases can coexist, because 7 C (A) < 7* (A), although the region, where 
both phases are stable is very narrow, since 7 C (A) and 7* (A) are numerically closed to each 
other. The phase transition is of the first order and takes place at some 70(A) between 
7c(A) and 7* (A) (fig. |2|). The line of phase transitions is defined as that line where the 



free energies of the both phases are equal to each other. Substituting po(9) into eq. (|2.4|) 
one can find the free energy per unit volume, to leading order in the large N limit: 

F J 0, in confining phase, (2 22) 



VN 2 \ ^E — \\c\ — ^g, in deconfining phase. 
The equations determining the critical point can be solved numerically, the result is: 

70(A) =4.219 A. (2.23) 
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3 Strong coupling lattice model in any dimension 



The model fll.31 ) is an unitary matrix analog of the Kazakov-Migdal model [20], which 



can be treated at large iV by the saddle point methods pi|-|2"4"|. Similar methods have 
been applied to the model ( |1.31| ) with A = PJ-0. Here we generalize the consideration 
of [|I[] to the case of A > 0. 

To begin, we choose the gauge in which the Polyakov loops are diagonal: g(x) = 
diag (^e tak ( x ^. Then, integrating over the gauge fields JJ\ we obtain an effective action for 
the eigenvalues ak(x): 



S, 



E 



+ 5> 

k<j 



. 2 0£ k {X) ~ 0£j(X) 



sin 



D 



+ o E lnl(a(x),a(x + n) [7) 



li=-D 



(3.1) 



where the second term comes from the Faddeev-Popov determinant and I(a, a' |r) is a 
one-link integral: 

I{a,a'\r)= J DU K{e ia ,U e ia 'U ] |r), (3.2) 

and we have used the invariance properties of the heat kernel. The link integral is calcu- 
lable for any iV [EH] and can be represented in the following form: 



I(a, a' |r) = const 



2tt 



2ttN 



, J(a) = ]^[ sin 



— a. 



(3.3) 



i<j 



J(a)J(a>) 

where d(z \r) is the Riemann theta function. 

In the large iV limit, since the effective action (|3.1|) is of order A^ 2 and there are A^ 



degrees of freedom, the statistical sum is dominated by the configuration which minimizes 
the action, i.e. the solution of the classical equation of motion. In terms of the eigenvalue 
density (|2.2|) the equation of motion reads: 



-2AIm(ci(a;)e^ e ) - p J* dO' p(9' \ x) cot 



9-9' 



6 



V — — 

J- D N 2 89 5p(9,x) 



lnl(a(x), a(x + //) [7). 



where 



c\{x) 



d9p(9,x)e 10 . 



(3.4) 



(3.5) 



The large A^ limit of the one-link integral ( p.2[ ) can be considered analogously to it's 



Hermitean-matrix counterpart along the lines of p6| , |24| . The method again is based 
on the correspondence between matrix quantum mechanics - cf. eqs. ( |1.29| ) and (|3.2|) 
(the integration in ( |3.2| ) acts as a projection on the singlet states) - and the quantum 
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mechanics of the free nonrelativistic fermions - note that the numerator in eq. fl3.3|) has a 
form of the Slater determinant. In the large N limit the fermions behave semiclassically 
and the problem reduces to the equations of motion of the collective field theory: 

da d , , 

* + »<«) = <>, (3-6) 

which should be solved on each link of the lattice with the following boundary conditions: 

a(9,0;x,p) = p(9,x), (3.8) 

a(9,y,x, p) = p(9, x + p). (3.9) 

The right hand side of eq. ( |3.4|) can be expressed through the solution of fl3~6[) - ( |3.9|) as 
follows: 



-2AIm(ci(a;)e- <0 ) + (D-l)p J* d9' p(9' , x) cot 



9-9' 



D 

= J2 (s(9 1 0;x,p) - s(0, 7; a; -//,//)). (3.10) 

jx=-D 

The equations ( |3.5|) - ( p,10|) completely determine the large N dynamics of the model 
under consideration. The vacuum should be identified with the ^-independent solution. 
At D = 1 these equations have a continuum limit. One should introduce the lattice 
spacing a, recover the canonical dimensions of the couplings, i.e. rescale 7 — > cry and 
A — > aX, and take the limit a — > 0. As a result of this procedure, one obtains eqs. ( |2.6| ) - 
(£§ of sec. 00. 

The vacuum solution in the large 7 confining phase is dictated by the center group 
symmetry - po(9) = To find the spectrum of excitations we linearize the equations of 
motion around po(9): 

p(9, X ) = ^ + ^Y. c n (x)e me , < = c_ n , (3.11) 

v(9,t;x,p) = ^- + ^- an(r;x,p)e me , a* n = a_ n , (3.12) 
s(9, r; x, l i) = y £ Pn(r; x, p) e m \ (3* n = /?_„. (3. 13) 



The solution of eqs. ( |3.6|) , (p.7|) linearized in a n and (3 n reads: 

. . _1_ / \ nT / \ tit . . 

a n [r\ x, p) — a^(x,p) e 2 + a n (x, p)e (3.14) 



. . % " 1 . . riT . . TIT , . 

p n (r; x, p) = - [a^(x, p) e 2 -a n (x,p)e 2 . (3.15) 
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Substituting these equations in the boundary conditions Q3.8Q , (|3.9| ) we express a n and 
P n in terms of c n : 



e T 2 Cn(x) - c n (x + fl) 



2 sinh ^ 



(3.16) 



The equality ( P-10Q then gives an equation for the Fourier coefficients of the eigenvalue 
density: 



D 



cosh - f — 1 + 5 n iX\ sinh ^ 



c n {x) + c n (x- n) 



which leads to the following mass spectrum: 



Ml = 2D (cosh ^ - l) - 2(D - 1 + <5 n i_DA) sinh ^. 



0, 

(3.17) 
(3.18) 



The strong coupling solution becomes unstable when M\ turns to zero. Form eq. ( 3.18|) 
one obtains: 

^ + y/p-l) 2 + AD[2(/J-l) + A/J] 
7c (A) = 21n J _ — . (3.19) 

We have not been able to obtain an exact solution in the weak coupling phase, but 
an approximate one, valid in the limit 7 — ► 0, can be found. At small 7 one can expand 



g(x) = e'v^M w 1 + z*/7$(a;), and the partition function ( 1.3 lj ) reduces to that of the 



Kazakov-Migdal model with the quadratic potential. From the point of view of the large 
iV equations of motion this approximation corresponds to the expansion of the left hand 
side of eq. (]3. ldf ) in the powers of 9. Actually, since at small 7 the eigenvalues are peaked 
about zero with the width of the distribution of order 7, we can expand the second term 
in eq. (|3.10| ), retaining only the contributions of order 7 _1 and 7 : 



p I d9' p{9',x)cot e -— - = 2p [ d9' pi - 9 ' ,X 



1 



J^F-ls 6 - ^ 



(3.20) 



The first term in ( |3. 10| ) is equal, with the same accuracy, to 2X9, as c\ = 1 + 0{p/ 2 ). In 
this approximation, we obtain the equations of motion for the Kazakov-Migdal model 
with the quadratic potential, as expected, and effective mass is equal to 



m 



off 



2D 



2\-~{D-l) 



7- 



(3.21) 



The vacuum solution of the latter model is known [22 



Po(0) 



71 



-fx 2 9 2 , 



myD - 1) + - 4(2D - 1 

(2D - 1)7 



(3.22) 
(3.23) 
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There are two reasons for which the weak coupling solution can terminate. First, it 
may become unstable due to the appearance of the massless excitation in the spectrum 
of fluctuations around it. For the solution ( 3.22[) , ( |3.23| ) of the Kazakov-Migdal model it 



happens at m 2 cS = 2\J2D — 1 ||23|| ; for smaller values of this solution does not exist 
since \i in ( 3.231) becomes complex. However, the analyses of the A = case shows || 



that the weak coupling phase terminates before it reaches the point of instability, because 
the eigenvalue distribution hits ir and begins to overlap with itself for smaller 7. More 
precise criterion for a critical point is the condition that the width of the support of a (9, r) 
exceeds 27r at some intermediate r between and 7 0. For the solution (|3.22|) , ( |3.23|) , 
the last criterion is satisfied when [[3]: 



The solution ( |3.22| ), fl3.23| ) being an approximate one, eq. ( |3.24| ) gives only an estimate 



for the critical coupling. In fact, it should give an upper bound. The arguments for 
this are the following: All higher terms in the expansion of the cotangent in (|3.20|) are 



negative, or, equivalently, the effective potential generated due to nonlinearity of the field 
g(x), is upside-down. Thus the neglected corrections can only strengthen the instability. 
These arguments are rigorous for A = 0. When A > 0, the first term in ( |3.10| ) which 
is proportional to A has an alternating Taylor expansion in powers of 9 when c\ is real. 
However, the first correction is also negative and we expect that (|3.24| ) gives an upper 
bound for the actual value of 7* (A) for all A. 

The consideration above gives the following picture of the phase transition (fig. |3[) . In 
all cases the confining and deconfining phases can coexist and are separated by a first- 
order phase transition. The model also undergoes a third-order large- N phase transition 
in the deconfining phase. Actually, the line on which the weak coupling solution termi- 
nates, 7 = 7* (A), obtained from the approximate equation (|3.24|) , crosses the line of the 



instability of the strong coupling phase 7 = 7 C (A) (fig. |3|). As ( |3.24j ) should give an upper 
bound for 7* (A), this crossing is not an artifact of the approximation done. So the weak 
coupling solution terminates before the transition to the strong coupling region and there 
exist two deconfining phases. The correlation length does not turn to zero at the critical 
line separating these phases. The phase transition is connected with the large- N critical 
behaviour of the link integral ( |3.2| ). Such third-order phase transitions are typical for 
unitary matrix models || IB, [T7J. 



The phase diagram of the D = 1 model is depicted on fig. |^a. One can verify that in 
the continuum limit 7 — ■> 0, A — > eq. ( |3.19| ) reduces to eq. ( |2.15| ) and eq. ( |3 . 24[ ) really 



gives an upper bound (7, (A) « ^ A = 48.70 A) for (ET2p 
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II 



continuum 

limit 21n(2D-l) 



A 



l/D 



a 



Figure 3: The phase diagram of the strong coupling lattice model: a - D — 1, b - D > 1. 
I - strong coupling (confining) phase, II, III - weak coupling (deconfining) phases; 1 - 
line on which the weak coupling phase terminates: 7 = 7* (A), 2 - line of the first-order 
phase transition, 3 - line of the instability of the strong coupling phase: 7 = 7 C (A), 4 - 
line of the third-order large N phase transition: 7 = 7* (A) 



4 Results and discussion 



We have considered the thermodynamics of the gas of quarks in adjoint representation of 
the gauge group interacting via non-Abelian electric forces. In 1+1 dimension the model 
that we have considered is adjoint QCD in the limit where quarks are heavy. The fugacity 
parameter is proportional to the Boltzmann weight of the classical particle of mass m - 
A oc e" m / T . The pre-exponential factor is not determined by classical theory, but it can 
be computed from a loop diagram: 



For the classical thermodynamics to be applicable, particle mass should be much larger 
than temperature - m » T. We also assume that m? » e 2 N and pair production 
is suppressed. The results of Sec. |2] indicate that there exists a region of parameters in 
which both of the conditions are satisfied and the model undergoes the deconfining phase 
transition. It is of the first order with critical line given approximately by the equation 



4.2^6-^-, (4.26) 
2T C V 2vr v ! 

or, explicitly, 

771 

T c = s 7 jT' (4-27) 

ln^ + 0(lnln^) 
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The value of T c determines, by standard arguments, the asymptotics of the density of 
states in 1+1 dimensional QCD with heavy adjoint matter - p(E) oc e E ^ Tc for large 
E. The spectrum and other properties of 1+1-dimensional QCD with adjoint matter 
fields have been investigated recently |27| - [51]. As emphasized by Kutasov | [2"7|] , this 



model shares some features with string theories in that it exhibits an infinite number of 
asymptotically linear Regge trajectories and its density states increases exponentially at 
high energy. It should therefore exhibit a Hagedorn temperature, which could either be 
an upper limiting temperature, or the critical temperature of a phase transition. Kogan 
and Zhitnitsky |30| outlined features which the spectrum would have to possess in order 
that the behavior is a phase transition. 

In the multidimensional case, we consider the limit of the lattice theory in which the 
coupling constant and the temperature measured in lattice units are large. In this limit 
the magnetic interactions can be neglected and the model also appears to be explicitly 
solvable. The deconfining phase transition for D > 1 takes place even in pure gluody- 
namics. For the model under consideration, this corresponds to A = 0, the case which 
was considered previously [0 [§]. The new feature of the model with adjoint matter 
is appearance of the additional third-order phase transition of Gross-Witten type for 
sufficiently large value of the fugacity parameter. 
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